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Abstract. Expressions for Sudakov form factors for heavy quarks are presented. They are used to construct 
resummed jet rates in e + e~ annihilation. Predictions are given for production of bottom quarks at LEP 
and top quarks at the Linear Collider. 
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1 Introduction 



momentum k± smaller than the mass, the collinear diver- 
gences are regularized by the quark mass. In this region 
mass effects are enhanced by large logarithms ln(mb/fcj_), 
increasing the significance of the measurement. Indeed, 
this leads to a multiscale problem since in this kinematical 
region also large logarithms \n(i/s/k±) appear such that 
both logarithms need to be resummed simultaneously. A 
solution to a somewhat similar two-scale problem, namely 
for the average sub-jet multiplicities in two- and three-jet 
events in e + e _ annihilation was given in We report 
here on the resummation of such logarithms in the fcj_- 
like jet algorithms ^2] and provide some predictions for 
heavy quark production. A preliminary comparison with 
next-to-leading order calculations of the three-jet rate [B] 
IT31ll4lll5| is presented. 



The formation of jets is the most prominent feature of 
perturbative QCD in e + e~ annihilation into hadrons. Jets 
can be visualized as large portions of hadronic energy or, 
equivalently, as a set of hadrons confined to an angular 
region in the detector. In the past, this qualitative defi- 
nition was replaced by quantitatively precise schemes to 
define and measure jets, such as the cone algorithms of the 
Weinberg-Sterman pQ type or clustering algorithms, e.g. 
the Jade or the Durham scheme (k± scheme) 0]. A 
refinement of the latter one is provided by the Cambridge 
algorithm Equipped with a precise jet definition the 
determination of jet production cross sections and their 
intrinsic properties is one of the traditional tools to in- 
vestigate the structure of the strong interaction and to 
deduce its fundamental parameters. In the past decade, 

precision measurements, especially in e + e~ annihilation, 2 Jet rates for heavy quarks 

have established both the gauge group structure underly- 
ing QCD and the running of its coupling constant a s over 
a wide range of scales. In a similar way, also the quark 
masses should vary with the scale. 

A typical strategy to determine the mass of, say, the 
bottom-quark at the centre-of-mass (cm.) energy of the 
collider is to compare the ratio of three-jet production 
cross sections for heavy and light quarks [6ll7ll8ll9lll()| . At 
jet resolution scales below the mass of the quark, i.e. for 
gluons emitted by the quark with a relative transverse 
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A clustering according to the relative transverse momenta 
has a number of properties that minimize the effect of 
hadronization corrections and allow an exponentiation of 
leading (LL) and next-to-leading logarithms (NLL) fj] 
stemming from soft and collinear emission of secondary 
partons. Jet rates in fc_L algorithms can be expressed, up 
to NLL accuracy, via integrated splitting functions and 
Sudakov form factors fj]. For a better description of the 
jet properties, however, the matching with fixed order cal- 
culations is mandatory. Such a matching procedure was 
first defined for event shapes in |16| . Later applications 
include the matching of fixed-order and resummed expres- 
sions for the four-jet rate in e + e~ annihilation into mass- 
less quarks |17U18| . A similar scheme for the matching of 
tree-level matrix elements with resummed expressions in 
the framework of Monte Carlo event generators for e + e~ 
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processes was suggested in ^5] and extended to general 
collision types in [201 - 

We shall recall here the results obtained in for 
heavy quark production in e + e~ annihilation. In the quasi- 
collinear limit (211122) . the squared amplitude at tree-level 
fulfils a factorization formula, where the splitting func- 
tions P a b for the branching processes a — > b + c, with at 
least one of the partons being a heavy quark, are given in 
D = 4 — 2e dimensions by 



p qq( z i<i) = C 'f 

P 9q( z ><1) = T R 



1 



e(l-z)- 



2z(l ~ z)m 2 



1 — z ' q 2 + (1 — z) 2 m 2 

22(1-2) 2z(l-z)m 2 



1 - 



1 



(1 - e)(q 2 + m 2 ) 



, (1) 



where z is the usual energy fraction of the branching, and 
q 2 is the space-like transverse momentum. As expected, 
these splitting functions match the massless splitting func- 
tions in the limit m — > for q 2 fixed. Furthermore, for 
q <C (1 — 2) m the splitting function Pqq is not any more 
enhanced at z — > 1, which is the known "dead cone" [2"3*| 
effect. The splitting function 
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+ 2(1-2) 



(2) 



obviously does not get mass corrections at the lowest or- 
der. 

Branching probabilities are defined through [15] 



1-9/O 

r Q (Q,q,m)= J dz 
q/Q 
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3 q 2 + m 2 

p g {Q,q)= J dzP gg (z)=2C A (lnQ-^j ,(3) 

9/Q 

and the Sudakov form factors, which yield the probability 
for a parton experiencing no emission of a secondary par- 
ton between transverse momentum scales Q down to Qo, 
read 



Aq(Q,Qo) = exp 



Q 



dq a s (q) 
q 7T 



r Q (Q,q) 



Ag(Q,Qo) = cxp 



On 



dq a s (q) 
q 7T 



(r g (Q,q) + r f (Q,q)) 



A f (Q,Q ) = [Aq(Q,Q )] /A g (Q, Q Q 



(4) 



where accounts for the number nj' h ^ of active 

light or heavy quarks. Jet rates in the k± schemes can 
be expressed by the former branching probabilities and 
Sudakov form factors. For the two- and the three-jet rates 



11 2 = {A Q (Q,Q )Y 



K 3 =2[A Q (Q,Q )Y 



(5) 
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dq a s (q) 
q 7r 



r Q (Q,q)A g (q,Q ) 



where Q is the cm. energy of the colliding e + e _ , and 
Qo = J/cutQ 2 plays the role of the jet resolution scale. 
Single-flavour jet rates in Eq. © are defined from the 
flavour of the primary vertex. 

In order to catch which kind of logarithmic corrections 
are resummed with these expressions it is illustrative to 
study the above formulae in the kinematical regime such 
that Q>m> Qo- Expanding in powers of a s , jet rates 
can formally be expressed as 



n2 
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k=n-2 



>(Q) 
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where the coefficients c 
L y = m(l/y cu t) and L n 
for the first order in a s are given by 



are polynomials of order I in 
ln(m 2 /<5Q). The coefficients 
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The coefficients for the second order in a s can be found 
in ^2], as well as the corresponding result for the four-jet 
rate. 

The impact of mass effects can be highlighted by two 
examples, namely by the effect of the bottom quark mass 
in e + e~ annihilation at the i?-pole (Fig. ^ up), and by 
the effect of the top quark mass at a potential Linear Col- 
lider operating in the TeV region (Fig.HJ down). In Fig.^ 
leading order (LO) and next-to-leading order (NLO) pre- 
dictions for three-jet rates are compared with the NLL 
result as obtained by numerical integration from Eq. 
For a comparison with the dead cone approximation see 
also |12| . While in the case of bottom quarks at LEP1 en- 
ergies the overall effect of the quark mass is at the few-per- 
cent level, this effect becomes tremendous for top quarks 
at the Linear Collider. Fixed order predictions for 6-quark 
production clearly fail at very low values of y cu t , by giving 
unphysical values for the jet rate, while the NLL predic- 
tions keep physical and have the correct shape. The latter 
is an indication of the necessity for performing such kind 
of resummations. Fixed order predictions work well for 
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Fig. 1. Effect of a &-mass of 5 GeV on the single- flavour 
three-jet rate at LEP1 energies as a function of the jet res- 
olution parameter in the Cambridge algorithm (upper plot). 
Effect of a t-mass of 175 GeV on the single-flavour three-jet 
rate as a function of the jet resolution parameter, at a poten- 
tial Linear Collider operating at cm. energies of 1 TeV (lower 
plot). Predictions are provided at LO, NLO and NLL accuracy. 

top production at the Linear Collider, a consequence of 
the strong cancellation of leading logarithmic corrections, 
and are fully compatible with our NLL result. 



3 Conclusions 

Sudakov form factors involving heavy quarks have been 
employed to estimate the size of mass effects in jet rates 
in e + e~ annihilation into hadrons. These effects are size- 
able and therefore observable in the experimentally rel- 



evant region. A preliminary comparison with fixed order 
results have been presented, and showed good agreement. 
Matching between fixed-order calculations and resummed 
results is in progress |24j . 



Acknowledgements 

We would like to thank Z. Trocsanyi for his kind invita- 
tion to present these results at the HEP2003 Europhysics 
Conference in Aachen. G.R. acknowledges partial support 
from Generalitat Valenciana under grant CTIDIB / 2002 /24 
and MCyT under grants FPA-2001-3031 and BFM 2002- 
00568. 



References 

1. G. Sterman and S. Weinberg, Phys. Rev. Lett. 39 (1977) 
1436. 

2. W. Bartel et al. [JADE Collaboration], Z. Phys. C 33 
(1986) 23. 

3. S. Bethke et al. [JADE Collaboration], Phys. Lett. B 213 
(1988) 235. 

4. S. Catani, Y. L. Dokshitzer, M. Olsson, G. Turnock and 
B. R. Webber, Phys. Lett. B 269 (1991) 432. 

5. Y. L. Dokshitzer, G. D. Leder, S. Moretti and B. R. Web- 
ber, JHEP 9708 (1997) 001 hep-ph/9707323 . 

6. G. Rodrigo, A. Santamaria and M. S. Bilenky, Phys. Rev. 
Lett. 79 (1997) 193 hep-ph/9703358 . 

7. P. Abreu et al. [DELPHI Collaboration], Phys. Lett. B 418 
(1998) 430. 

8. A. Brandenburg, P. N. Burrows, D. Muller, N. Oishi and 
P. Uwer, Phys. Lett. B 468 (1999) 168 hep-ph/9905495 . 

9. R. Barate et al. [ALEPH Collaboration], Eur. Phys. J. C 
18 (2000) 1 [hep-ex/00080i"3l . 

10. G. Abbiendi et al. [OPAL C ollaboration], Eur. Phys. J. C 
21 (2001) 411 hep-ex/0105046 . 

11. S. Catani, B. R. Webber, Y. L. Dokshitzer and F. Fiorani, 
Nucl. Phys. B 383 (1992) 419. 

12. F. Krauss a nd G. Rodrigo, |hep-ph/0303038l 
hep-ph/0309326 

13. G. Rodrigo, Nuc l. Phys. Proc. Suppl. 54A (1997) 60 
hep-ph/9609213 . 

14. M. S. Bilenky, S. Cabrera, J. Fuster, S. Marti, G. Ro- 
drigo and A. Santamaria, Phys. Rev. D 60 (1999) 114006 
hep-ph/9807489 . 

15. G. Rodrigo, M. S. Bilenky and A. Santamaria, Nucl. Phys. 
B 554 (1999) 257 hep-ph/990 5"276"l; Nucl. Phys. Proc. 
Suppl. 64 (1998) 380 hep-ph/9709313 . 

16. S. Catani, L. Trentadue, G. Turnock and B. R. Webber, 
Nucl. Phys. B 407 (1993) 3. 

17. L. J. Dixon and A. Signer, Phys. Rev. D 56 (1997) 4031 
hep-ph/9706285 . 

18. Z. Nagy and Z. Trocsanyi, Phys. Rev. D 59 (1999) 014020 
[Erratum-ibid. D 62 (2000) 099902] |hep-ph/9806317|. 

19. S. Catani, F. Krauss, R. Kulm and B. R. Webber, JHEP 
0111 (2001) 063 [hep-ph/010923"!] . 

20. F. Krauss, JHEP 0208 (2002) 015 hep-ph/0205283 . 

21. S. Catani, S. Dittmaier and Z. Trocsanyi, Phys. Lett. B 
500 (2001) 149 hep-ph/0011222 . 



4 German Rodrigo, Frank Krauss: Resummed jet rates for heavy quark production in e e annihilation 



22. S. Catani, S. Dittmaier, M. H. Seymour and Z. Trocsanyi, 
Nucl. Phys. B 627 (2002) 189 hep-ph/0201036 . 

23. Y. L. Dokshitzer, V. A. Khoze and S. I. Troian, J. Phys. 
G 17 (1991) 1602. 

24. F. Krauss, G. Rodrigo, A. Schalicke and P. Tortosa, in 
preparation. 



